Fine structure of the isoscalar giant quadrupole 
resonance in 40 Ca due to Landau damping? 

I. Usman a < b , Z. Buthelezi a , J. Carter b , G.R.J. Cooper c , R.W. Fearick d , 
S.V. F6rtsch a , H. Fujita a ' b , Y. Fujita , Y. Kalmykov f , 
P. von Neumann-Cosel f '*, R. Neveling a , P. Papakonstantinou f , A. Richter f ' s , 
R. Roth f , A. Shevchenko f , E. Sideras-Haddad b , F. D. Smit a 

o 

' a iThemba LABS, PO Box 722, Somerset West 1129, South Africa 

$_| , School of Physics, University of the Witwatersrand, Johannesburg 2050, South Africa 

, c School of Earth Sciences, University of the Witwatersrand, Johannesburg 2050, South 

Africa 

a Department of Physics, University of Cape Town, Rondebosch 7700, South Africa 
"Department of Physics, Osaka University, Toyonaka, Osaka 560-0043, Japan 
Hnstitut fur Kernphysik, Technische Universitdt Darmstadt, D-64289, Darmstadt, Germany 
SECT*, Villa Tambosi, 1-38123 Villazzano (Trento), Italy 

X 
=5 

a 



(N 
On 
O 

cn 
o 



Abstract 

The fragmentation of the Isoscalar Giant Quadrupole Resonance (ISGQR) in 
40 Ca has been investigated in high energy- resolution experiments using proton 
inelastic scattering at E p = 200 MeV. Fine structure is observed in the region 
of the ISGQR and its characteristic energy scales are extracted from the ex- 
perimental data by means of a wavelet analysis. The experimental scales are 
well described by Random Phase Approximation (RPA) and second-RPA cal- 
culations with an effective interaction derived from a realistic nucleon-nucleon 
interaction by the Unitary Correlation Operator Method (UCOM). In these re- 
sults characteristic scales are already present at the mean-field level pointing to 
their origination in Landau damping, in contrast to the findings in heavier nu- 
clei and also to SRPA calculations for 40 Ca based on phenomenological effective 
interactions, where fine structure is explained by the coupling to two-particle 
, two-hole (2p-2h) states. 
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1. Introduction 

High energy-resolution experiments on nuclei excited in the region of giant 
resonances present a unique way to extract information about the dominant 
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processes of the decay. This information is imprinted into the fragmentation of 
giant resonances visible as fine structure in excitation energy spectra. Recent 
work has established fine structure as a global phenomenon in medium-mass to 
heavy nuclei for the case of the Isoscalar Giant Quadrupole Resonance (ISGQR) 
[1] . Fine structure has also been demonstrated for a variety of other modes like 
the spin-isospin flip Gamow-Teller (GT) mode 0, the Isovector Giant Dipole 
Resonance (IVGDR) and the magnetic quadrupole resonance @. 

A variety of theoretical approaches has been put forward to understand the 
fine structure including doorway-state analysis 0, a local scaling dimension 
model Q , an entropy index method [§] and wavelet analysis [HI El ■ A com- 
parison for representative cases indicates that a wavelet analysis is a particularly 
promising tool since it provides simultaneously a quantitative measure of 
the fine structure and information on the localization in the excitation energy 
spectrum. Characteristic scales can be extracted from the power spectra of 
wavelet transforms, which allow a direct comparison between experiment and 
theory. A systematic study of the ISGQR |l[ shows that the observed scales in 
medium-mass to heavy nuclei originate from a collective damping mechanism 
induced by the coupling of elementary one-particle one-hole (lp-lh) states to 
low-lying surface vibrations 13j. 

Recently the high energy-resolution measurements of the ISGQR have been 
extended 14| to the low-mass region 12 < A < 40. These experiments were 
motivated by the question whether collective damping remains the most impor- 
tant decay mechanism in light nuclei or whether other contributions like direct 



decay or Landau damping become more important [15|. The present work fo- 
cuses on the case of 40 Ca where considerable fragmentation of E2 strength has 
been observed in electron [§], proton [l6, 17 1 and alpha scattering [l|| studies, 
albeit measured with varying resolution. The choice of a doubly magic nucleus 
allows a comparison to calculations within the framework of RPA and SRPA 
and thus to distinguish between the role of lp-lh and 2p-2h excitations. Since 
in the calculations a modern realistic interaction derived b y t he Unitary Cor- 
relation Operator Method (UCOM) described in Refs. 1^-21| has been used, 
the experimental data serve at the same time as a test bench for UCOM. This 
interaction is one of a family of realistic interactions and has been used with 
quite some success in the description of basic nuclear quantities such as masses 



22j and the gross features of giant resonances, in particular the giant dipole 



and giant quadrupole ones on the SRPA level |23j, \2A \ . 



2. Experiment 

The experiments were carried out with a 200 MeV proton beam produced 
by the Separated Sector Cyclotron (SSC) of the iThemba Laboratory for Accel- 
erator Based Sciences (iThemba LABS), South Africa. The protons scattered 
inelastically on a 40 Ca target with an areal density of 3.0 mg/cm 2 were momen- 
tum analyzed with a K600 magnetic spectrometer after scattering. Dispersion 
matching techniques were used in order to exploit the high energy-resolution 
capability of the spectrometer. Energy resolutions AE — 35 — 40 keV Full 
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Width at Half Maximum (FWHM) were achieved. The scattering angles were 
selected to be below, at and above the maximum of the cross section for AL = 2 
transitions into the ISGQR. Therefore, measurements were performed at three 
different scattering angles, #Lab = 7°, 11° and 15°, as illustrated in Fig.[TJ for ex- 
citation energies between 6 and 30 MeV. Details of the data analysis procedures 
are described elsewhere 
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Figure 1: Excitation energy spectra for 40 Ca at scattering angles of $Lab = 7°, 11° and 15°. 
Note that the ISGQR is most strongly excited at $Lab = 11°. The drop of counts in the 
spectra above 27 MeV is due to a gradual loss of angular acceptance towards the end of the 
spectrometer focal plane. 



The excitation energy spectrum in Fig. [T]taken at #Lab = 11°, which corre- 
sponds to the maximum cross section for the ISGQR, reveals a broad resonance 
at a mean energy £ x « 18 MeV with strength distributed between approxi- 
mately 12 and 22 MeV. Intermediate structure is visible through peaks around 
12 , 14, 16, 17 and 18 MeV consistent with previous experimental work 
Il8j |. Pronounced fine structure is observed up to about 20 MeV. At the larger 
scattering angle of #Lab = 15°, the overall structure is reduced in size but dis- 
tinctive features of the fine structure persist. At the smaller scattering angle 
$Lab = 7°, the fine structure changes somewhat indicating the presence of other 
multipoles. 

The dominant quadrupolc nature of the bump on top of a smooth quasifree 
background has been demonstrated by multipole decompositions of angular 
distributions measured in alpha [26[ and proton [l6| scattering. Properties of 
the ISGQR have also been investigated with electromagnetic probes [a, |27| and 
capture experiments [28] . The results are consistent in that the energy weighted 
sum rule for isoscalar E2 strength is almost exhausted in the excitation range 
of interest. Strong transitions at energies E K < 12 MeV consist of a mixture of 
different multipoles [29[ , while the peaks in the excitation region around 14 MeV 
were shown to have quadrupole character [301 ] . 
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The multipole analysis of proton scattering data in Ref. [16[ indicates that 
for the kinematics of the present experiment at the optimum angle for AL = 2 
transfer, contributions from the overlapping isoscalar giant monopole resonance 
(ISGMR) and low-energy AL = 3,4 strength are at the level of a few percent 
at best. The IVGDR is strongly populated at forward scattering angles by rel- 
ativistic Coulomb excitation. However, the angular distribution falls off steeply 
such that the contribution is already small (but still sizable) at 6*Lab = 7° and 
negligible at 11°. Furthermore, an analysis of the spectra along the lines de- 
scribed in Ref. Q with the additional constraint of a constant level density 
independent of the kinematics indicates only limited contributions to the cross 
sections in the ISGQR region from quasifree processes and/or other multipolcs. 
Thus, it seems justifed to assume that the fine structure of the cross sections at 
@Lab — 11° arises solely from the ISGQR. 

Another way to verify this conjecture is to use a cross correlation analysis 
as described e.g. in Ref. [3l[. For that purpose, any smooth background in 
the spectra was removed by a discrete wavelet analysis @, [l2[ • In a next step 
gross structures were taken out by folding with a Gaussian of appropriate width. 
The ratio between these two spectra - the so-called stationary spectrum - then 



contains local fluctuations around the mean only 32]. The normalized cross 
correlation functions Ce 1 ,e 2 between stationary spectra including the 7° angle 
are quite similar: C70 n° = 0.26(8) and CV°,i5° = 0.21(7). In contrast, a much 
larger value Cno il5 o = 0.50(8) is found for the 11°, 15° angle pair indicating a 
high degree of correlation, since the value is above an expression derived for the 
angular cross-correlation function of cross-section fluctuations [33| . We conclude 
that the fine structure observed in the 11° spectrum is indeed due to the ISGQR 
and contributions from other modes must be small. 



3. Extraction of scales 



Wavelet analysis techniques |34| have been utilized for the extraction of 
characteristic energy scales from the region of the ISGQR excitation [l[. Con- 
tinuous Wavelet Transforms (CWT) using both the real and complex Morlet 
mother wavelets have been successfully applied previously [l2| ■ Application of 
the CWT to the present data in the energy region of the ISGQR is shown in 
Fig.O By plotting the real part of the complex coefficients in a two-dimensional 
distribution of energy scales vs. excitation energy, the positions of the structures 
within the original energy spectrum can be identified. Maxima of the wavelet 
coefficients at certain scale values over the energy region of ISGQR (or parts of 
it) indicate characteristic scales. This is illustrated in the middle part of Fig. [2l 
where the plotted wavelet energy scale was restricted to 3.0 MeV in order to 
reveal scales related to intermediate structure. The bottom part of Fig. [2] shows 
the same results but with an expanded wavelet scale to emphasize the presence 
of characteristic scales in the fine structure as well. 

In order to obtain a quantitative measure for the characteristic scales, the 
absolute values of the complex coefficients are projected onto the wavelet scale 
axis. The resulting power spectrum obtained is shown in r.h.s. of Fig. E^a) 
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Pi gure 2: Excitation energy spectrum and wavelet transform for ^Oa.(j?,p') at $Lab — 
using the complex Morlet mother wavelet for scales up to 3.0 MeV. The bottom part shows 
the presence of characteristic scales in the fine structure below 800 keV. Characteristic scales 
identified by maxima in the power spectrum are indicated by dashed lines. 



and maxima representing characteristic scales are indicated by arrows. Three 
peaks at scale values 150, 240 and 460 keV arising from the fine structure can 
be identified as well as three scales at 1.05, 2.0 and 3.9 MeV corresponding to 
intermediate structure. 



4. Model calculations and discussion 



To the best of our knowledge, all existing RPA calculations of the ISGQR 
in 40 Ca (see Refs. [35l - l37j l for some recent examples) except the one discussed 
below and shown in Fig.[3jb) on the l.h.s. are characterized by a dominant peak, 
which cannot yield a scale in the wavelet power spectrum (cf. Fig. 13 in Ref. [H). 
We have verified this result using Skyrme and Gogny interactions. Thus, the 
experimentally observed fragmentation is usually attributed to mixing of the 
lp-lh states with more complex configurations, including collective phonons 

As such, here we perform large-scale RPA and SRPA calculations with a re- 
alistic interaction, Vtjcom- The main differences of the present to other SRPA 
approaches are: (i) The two-body, finite-range interaction is the sole input, de- 
scribing both the ground state and the residual couplings. In particular, we use 
Vucom, a- potential derived from the Argonne V18 interaction by renormalizing 
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Figure 3: L.h.s.:Experimental spectrum of the 40 Ca(p,p') reaction at the maximum cross 
section of the ISGQR (a) and theoretical isoscalar quadrupole strength distributions in 40 Ca 
from RPA (b), full SRPA (c), SRPA including 2p-2h transitions below 200 MeV only (d) and 
SRPA including 2p-2h transitions above 200 MeV only (e), based on the UCOM interaction. 
R.h.s.: Corresponding power spectra from the CWT. Characteristic scales are indicated by 
arrows. 



it within the Unitary Correlation Operator Method |20|, |2l|. (ii) No trunca- 
tion is imposed on the model space describing the excited states, other than 
that due to the exhaustion of the single-particle space. This means that the 
Hartree-Fock (HF) equations are solved within a selected single-particle space 
and all lp-lh and 2p-2h configurations available within that space are taken into 
account when building the SRPA matrix. The single-particle basis used here 
consists of harmonic-oscillator wave functions with orbital angular momentum 
£ < 6 and node quantum number n < 6, providing reasonable convergence. The 
number of lp-lh and 2p-2h configurations is about 10 6 . Further details on the 



SRPA approach can be found in Refs. [23l. |24|. 



The isoscalar quadrupole strength distributions and power spectra resulting 
from both RPA and SRPA are shown in Figs. [3th) and (c), respectively. While 
the RPA results using Vucom sytematically overestimate ISGQR energies [HI, 
a striking result is that the ISGQR strength in 40 Ca already appears fragmented 
at the RPA level. The same interaction produces one dominant ISGQR peak for 
heavier nuclei like 90 Zr and 208 Pb 43j. We return to this point later. In SRPA 
the resonance is shifted strongly to lower energies compared to RPA. Indeed, 
since Vucom is a bare interaction with respect to long-range correlations, the 
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SRPA is employed not only to supplement RPA with a damping mechanism, 
but also to effectively dress the Lp-lh configurations, affecting their energy (for 



a relevant discussion see Refs. [23|, |24|). Also observed is that the resonance 



strength and the profile of the strength distribution are altered, but, as is usually 
the case with this type of calculation 2J, |44j , not dramatically. The SRPA with 
the present interaction reproduces the global properties of the ISGQR well (the 
same holds for the IVGDR), considering that the input has not been fitted to 
this end [2J| . In this work it is used for the first time to analyze the fine structure 
of the ISGQR. 

The power spectra analysis reveals that the intermediate and large scales are 
comparable in RPA and SRPA, but in SRPA small scales are somewhat more 
important in the power spectrum. Characteristic scales indicated by arrows in 
Fig.[3]show a good correspondence between the theoretical results (for both RPA 
and SRPA) and the experimental data except for the smallest scale where the 
relation to the calculations is unclear, and the largest one, which is suppressed 
in the calculations. 

The question arises as to the origin of the various scales. The RPA model ac- 
counts only for Landau damping. Therefore, our RPA results, Fig. [3Kb), suggest 
that this mechanism plays an important role in the case of Ca (but not in heav- 
ier nuclei), contrary to predictions using phenomenological effective interactions 
and density functionals, as already mentioned. Regarding the \pti~ 1 ) structure 
of the dominant peaks, we find that the leading components as well as their 
spacings are quite similar in RPA and SRPA. For example, the two strongest 
peaks at the higher-energy part of the resonance (RPA: 26.81, 28.75 MeV; SRPA: 
17.17, 19.09 MeV) in both cases contain large \0g, 0d _1 ) components, while the 
lowest prominent peak (RPA: 21.59 MeV; SRPA: 14.88 MeV) in both cases is 
rather dominated by \ld, (Odor Is) -1 ) configurations. 

Apparently, Landau damping accounts for most of the scales. This result 
must be a consequence of the properties of the Vucom, which differs in many 
respects from phenomenological effective interactions. Firstly, it retains all the 
complexities of the realistic NN interaction and has not been fitted using HF 
and RPA calculations (therefore it yields, e.g., underbound nuclei in HF and too 
high GQR energies in RPA calculations). We found strong Landau fragmenta- 
tion also with a different renormalizcd realistic interaction, namely the CD-Bonn 
potential transformed via the Similarity Renormalization Group (SRG), suggest- 
ing that we are not observing an artifact of the UCOM procedure. Semi-realistic 
effective interactions supplemented with a phenomenological three-body contact 
term (45j . leading to an improved description of nuclear radii (SB} and of the 
three major giant resonances on the RPA level still produce considerable 
fragmentation, by comparison. Regarding the tensor terms of Vucom which 
do not affect the HF ground state, we find only a weak influence on the RPA 
strength function, although in SRPA they account for roughly half of the en- 
ergy shift of the ISGQR between RPA and SRPA results. Secondly, the present 
interaction does not parameterize long-range correlations and corresponds to a 
very low nucleon effective mass m* /m in HF, as is evident from the very broad 
unperturbed spectra 21, 43|. This alone, however, seems not to be a critical 
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property. We find, for example, that for traditional Skyrme functional and 
using continuum RPA, a low m*/m value does not lead to fragmentation. Nei- 
ther does the coupling with the particle continuum produce sufficient spreading 
(only up to about 1 MeV FWHM). 

Next, we examine the role of 2p-2h configurations in the vicinity of the 
resonance and of high-lying ones (the present model space contains 2p-2h con- 
figurations with energies up to about 600 MeV). In Figs. [3][d) and (e), SRPA 
results are shown using only 2p-2h configurations with energies below or above 
200 MeV, respectively, the latter being roughly the median 2p-2h energy within 
the present model space. The configurations in the vicinity of the resonance 
alone tend to produce spreading, but the corresponding scales are somewhat 



too large and the energy of the resonance is also too high, as expected [24 1. 
Subsequent inclusion of the highest states compresses and shifts the spectrum 
to lower energies. Conversely, high-lying 2p-2h configurations alone tend to 
produce smaller scales. 

The smallest SRPA scale must reflect, to some extent, the structure of the 
density of states in the vicinity of the resonance and in particular the distri- 
bution of 2p-2h configurations acting as doorway states. Therefore, its precise 
value and prominence could be sensitive to the diagonal approximation (use of 
unperturbed 2p-2h configurations) employed in the present SRPA calculation. 
That approximation is expected to have no other significant influence on our 
results though 23|, [24 [ . 

Within the Quasiparticle Phonon Model (QPM) and the Extended Time 
Dependent Hartree-Fock (ETDHF) model, it has been possible to attribute 
scales observed in the ISGQR to the coupling to low- lying surface vibrations 
(collective mechanism) [l], In the spirit of Refs. [3, [13], we have attempted 
to solve the SRPA problem while excluding or selecting the strongest matrix 
elements coupling the Ip-Ih and the 2p-2h spaces. The probability distribution 
of the matrix elements here is similar to QPM results (strongly cusp-shaped 
around zero). No conclusion could be reached as to the effect of the few strongest 
or the many weak couplings. Small scales were present in all cases. We note that 
the present formalism contains only incoherent couplings to the 2p-2h space, 
although some collectivc-phonon effects may be present on a Tamm-Dancoff 
approximation level, to the extent that the diagonal approximation has been 
found to be rather good. In any case, it is not possible to isolate collective 
effects unambiguously, if at all. 

No clear signature is found in the calculatons for the smallest experimen- 
tal scale. Fine structure in the experimental spectra at these scales is most 
likely due to overlapping Compound nuclear states, i.e. Ericson fluctuations [47 [ 
damped by the experimental resolution 3l|. Finally, regarding the largest ex- 
perimentally extracted scale at about 4 MeV, we can only speculate that it is 
due to the smooth part of the resonance (its width), not reproduced by the cal- 
culations. The smooth part certainly contains contributions due to a coupling 
with the particle continuum, broadening all peaks. Escape widths in light nuclei 
up to 40 Ca have been shown to be significant [15j . Also higher-order configu- 
rations and correlations may play a role. We note that in some existing SRPA 
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and extended shell-model calculations empirical single-particle states are used, 
i.e., the input is already correlated which may be the reason why they appear 
more efficient in producing more or less a smooth spreading (39l . [iol |48| . This 
discussion, however, is beyond the scope of the present work. 



5. Conclusions 

Fine structure of the ISGQR in 40 Ca has been observed in high energy- 
resolution proton scattering experiments and analyzed in the framework of a 
wavelet approach. Comparisons with RPA and SRPA calculations using the 
UCOM interaction reveal that the scales characterizing the fine structure are 
largely present already on the mean- field level. This finding is in contrast to 
previous calculations of the ISGQR strength distribution in 40 Ca, which at- 
tribute the fragmentation to the coupling of Ip-Ih states to low-lying phonons 
38- 42j • At present we cannot exclude that the agreement of scales deduced 



from RPA and SRPA is - at least partly - fortuitous. However, tests using dif- 
ferent interactions indicate that the above result is possibly a consequence of the 
properties of the renormalized interaction used, which retains all the complexi- 
ties of the realistic NN interaction as compared to phenomenological mean-field 
interactions. The dressing of lp-lh states in a SRPA calculation is necessary 
to achieve a proper reproduction of the mean ISGQR energy and to improve 
details of the strength distribution. It seems unlikely that a full treatment of the 
particle continuum can lead to fragmentation. At the same time, we cannot rule 
out that coupling to collective phonons additionally plays a role in generating 
the experimentally observed scales. The present results underline the relevance 
of studies of the fine structure of giant resonances as a tool for an improved 
understanding of collective modes and their decay. 
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Note added in proof 

Recent work by Aiba et al. [HI], published after submission of our manuscript, 
circumstantiates our conclusions based on an extraction of the local scaling 
dimension of fine structure fluctuations in shell- model calculations of the 40 Ca 
GQR strength function. 
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